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Abstract 

A fractional generalization of variations is used to define a stability of non-integer 
order. Fractional variational derivatives are suggested to describe the properties of 
dynamical systems at fractional perturbations. We formulate stability with respect 
to motion changes at fractional changes of variables. Note that dynamical systems, 
which are unstable "in sense of Lyapunov", can be stable with respect to fractional 
variations. 

The theory of integrals and derivatives of non-integer order goes back to Leibniz, Liouville, 
Riemann, Grunwald, and Letnikov [JJ, [2] . Fractional analysis has found many applications in 
recent studies in mechanics and physics. The interest in fractional integrals and derivatives 
has been growing continuously during the last few years because of numerous applications. In 
a short period of time the list of applications has been expanding. For example, it includes 
the chaotic dynamics [3j 0], material sciences [5j EJ [7J, mechanics of fractal and complex 
media [H [9], quantum mechanics [TOj ITT] , physical kinetics [31 [121 021 ELU EES], long-range 
dissipation [H1[T7], non-Hamiltonian mechanics [T8| [T9]. long-range interaction [221 [23[ 124"]. 

In this preprint, we use the fractional generalization of variation derivatives that are 
suggested in [19]. Fractional integrals and derivatives are used for stability problems (see 
for example [251 I2H1 [271 I2H1 [29]). We consider the properties of dynamical systems with 
respect to fractional variations [19]. We formulate stability with respect to motion changes 
at fractional changes of variables (see pages 294-296). Note that dynamical systems, 
which are unstable "in sense of Lyapunov" , can be stable with respect to fractional variations. 

Let us consider dynamical system that is defined by the ordinary differential equations. 
Suppose that the motion of dynamical system is described by the equations 

^y k = F k (y), k = l,...,n. (1) 

Here yi, ...,y n be real variables that define the state of dynamical system. 

Let us consider the variation 8y k of variables y k . The unperturbed motion is satisfied to 
zero value of variation Sy k = 0. The variation describes that as function f(y) changes at 
changes of argument y. The first variation describes changes of function with respect to the 
first power of changes of y: 

5f(y) = Dlf(y)dy, (2) 
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where 

The second variation describes changes of function with respect to the second power of 
changes of y: 

5 2 f{y)=Dlf{y){dyf. (3) 

The variation 5 n of integer order n is defined by the derivative of integer order Dyf(y) = 
d n f/dy n , such that 

5 n f(y)=D n y f(y)(dy) n . 
Let us define the variation of fractional order [19] by the equation 



S a f = D a y f (5y) a , (4) 

where D° is a fractional derivative [H [2] with respect to y. 

The fractional variation of order a describes the function f(y) changes with respect to 
fractional power of variable y changes. The variation of fractional order is defined by the 
derivative of fractional order. 

Let us derive the equations for fractional variations 5 a y k - We consider the fractional 
variation of equation (CQ) in the form: 

* a 4l/fc = k = l,...,n. (5) 

at 

Using the definition of fractional variation (J3J), we have 

S a F h {y) = [D° l F k ](6y l ) a , fc=l,...,n. (6) 
From Eq. (Q, and the property of variation 

5 a ^-y k = ^8 a y k , (7) 
dt y dt y ' v ; 



where y k = yk(t, a), we obtain 

d 



dt 



6 a y k = iD«F k }(6 yi y } h = l,...,n. 



Note that in the left hand side of Eq. (151) , we have fractional variation of 5 a yk, and in the 
right hand side - fractional power of variation (5y k ) a . 

Let us consider the fractional variation of the variable y k . Using Eq. (jlj), we get 

S a y k = [D%y k ] (S Vl r, k = l,...,n. (9) 
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For the Riemann-Liouville fractional derivative, 

D«y h * 

ilk ^l. Therefore for simplification of our transformations, we use the fractional derivative 
fractional power of derivative 

D a yi y k = b kl Dlyu (10) 

where 5 k i is a Kronecker symbol. Substituting Eq. ffTUl) into Eq. (Q, we can express the 
fractional power of variations (Sy k ) a through the fractional variation 5 a y k : 

(6y k ) a = D a yk y k 5 a y k . (11) 

Substitution of Eq. ffTTT) into Eq. ([8]) gives 

d 



dt 



S a Vk = [D a yi yi] [D«F k ] 5<*y, (12) 



Here we mean the sum on the repeated index I from 1 to n. Equation ffl~2l) is equations for 
fractional variations. Let us denote x k the fractional variations 5 a y k : 

x k = 6 a y k = [D a yk y k ] (5y k ) a . (13) 

As a result, we obtain the differential equation for fractional variations 

^x k = a k i(a)x h (14) 

where 

o«(a) = [L>>] D a yi F k . (15) 

Using the matrix X 1 = (xi, x n ), and A a = \\a k i(a)\\, we can rewrite Eq. ( |T4"1) in the 
matrix form 

±X-AJC. (16) 

Equation (|T6l) is a linear differential equation. To define the stability with respect to frac- 
tional variations, we consider the characteristic equation 

Det[A a - XE) — (17) 

with respect to A. If the real part i?e[Afc] of all eigenvalues A for the matrix A a are negative, 
then the unperturbed motion is asymptotically stable with respect to fractional variations. 
If the real part -Re[Afc] of one of the eigenvalues A of the matrix A a is positive, then the 
unperturbed motion is unstable with respect to fractional variations. 
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A system is said to be stable with respect to fractional variations if for every e, there is 
a 5o such that: 

\\5 a y(t )\\ < 5 => \\5 a y(a,t)\\ < e Vi G R+. (18) 

The dynamical system is said to be asymptotically stable with respect to fractional variations 
5 a y(t, a) if as 

t^oo, \\S a y(t,a)\\^0. (19) 

The concept of stability with respect to fractional variations (see [20] pages 294-296) is 
wider than the usual Lyapunov or asymptotic stability. Fractional stability includes concept 
of "integer" stability as a special case (a = 1). A dynamical system, which is unstable 
with respect to first variation of states, can be stable with respect to fractional variation. 
Therefore fractional derivatives expand our possibility to research the properties of dynamical 
systems. 
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